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Quantitative predictions are presented of a depletion-induced torque and force acting on a single 
colloidal hard rod immersed in a solvent of hard spheres close to a planar hard wall. This torque and 
force, which are entirely of entropic origin, may play an important role for the key-lock principle, 
where a biological macromolecule (the key) is only functional in a particular orientation with respect 
' to a cavity (the lock). 
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C\) \ The depletion effect, i.e., the effect that smaller particles in a colloidal mixture are expelled from those regions where 
bigger particles are sufficiently close together, leads to an effective force between the bigger particles. In the case of 
hard-core interactions these effective forces are purely entropic in origin [QJ^]. For the relatively simple geometries of 
two large spheres or a sphere close to a planar or curved wall in a solvent of small hard spheres this depletion force 
has been studied in detail in theory simulation J|Q, and experiment J|,|§]. Depletion forces are crucial for the 

phase behavior of colloidal system; they can drive gas- liquid and liquid-solid phase separation, e.g., in colloid-polymer 
fiof and colloid-colloid mixtures (n|. It has been suggested that depletion effects play an important role in the 
biological "key-lock" mechanism, in which a nonspherical macromolecule (the key) fits into a cavity (the lock) and 
| forms a chemical bond, but only if the key has a particular orientation. A robust performance of this mechanism 
is possible if the biological environment is capable of passively transporting the key, correctly oriented, to the lock. 
This requires a chemically unspecific force to transport the center-of-mass of the macromolecule towards the cavity 
and a torque to orient it. In order to be able to study the relevance of the entropic contribution to this mechanism 
a quantitatively reliable tool for predicting such an entropic torque is needed. Using density functional theory and 
computer simulations, we show and predict quantitatively that depletion effects can generate both a force and a 
torque on a nonspherical particle in the vicinity of a planar wall thereby providing further evidence for the relevance 
of depletion phenomena in key-lock mechanisms. 
| The specific system we study is illustrated in Fig|j] and consists of (i) an infinite planar hard wall located at z = 0, 
£SJ . (ii) a hard spherocylinder with cylinder length L, diameter a, at distance z from the wall and forming an angle 
9 E [0, tt/2] with respect to the wall normal, and (iii) a solvent of hard spheres with diameter a s and number density 
p s in the bulk, i.e., far away from the wall. Since we are considering only hard interactions the temperature T of the 
system plays only the role of an energy scale via (3 — X/ksT. We are interested in the effective rod- wall potential 
W(z, 9) induced by the spheres. The force and the torque on the rod follow by differentiating —W(z, 9) with respect 
to z and 9, respectively. The bare rod- wall interaction restricts the center of mass of the rod to z to z > z m i n (9) = 
(a + L\ cos#|)/2. The system, and hence (3W(z, 9), is completely characterized by the aspect ratio L/a of the rod, 
' O ' the diameter ratio cr/a s , and the packing fraction rj — p s TTa^./d>. 

Following Asakura and Oosawa [Q, a first approximation of the effective potential W(z,9) can be obtained from 
the analysis of the volume that is excluded to the centers of the solvent spheres by the wall and the rod. The wall 
excludes the slab < z < cr s /2, and the rod the volume of a spherocylinder of length L and diameter a + a s . If 
these two contributions to the excluded volume overlap, i.e., for z < z m i n + a s (see Fig. Q) the volume accessible to 
the centers of the spheres and hence the entropy of the solvent increases. This gain of entropy translates directly 
into an effective, purely attractive force acting on the center of the rod. An important ingredient of this so-called 
Asakura-Oosawa approximation (AOa) is that the density profile of the spheres is taken to be constant and equal to 
the bulk density p s , i.e., all correlations in the solvent are entirely neglected. The AOa therefore underestimates the 
strength and the range of the force acting on the rod. 

Here we go beyond the AOa taking the sphere-sphere correlations into account. We recall that the depletion 
potential equals the grand-potential change of the (inhomogeneous) fluid of spheres in contact with the wall upon 
moving the rod from z — > oo to a finite value of z, at a given 9. Within density functional theory (DFT), with 
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Tex [ps , Pr] the excess (over ideal) free energy functional of a mixture of hard spheres and rods, one can write [E| 

0W(z, 0) = lim [c r 1] (z-kx>,8)- 4 1} (z, 9)} , (1) 

where c£ (z, 0) = — 5/3^-" e a:|/ > s! Pr]/3pr{z-i is the direct one-body correlation function of the rods |D| where p r (z,9) 
is the number density of rods for a given orientation 9. The functional we use is based on Rosenfeld's fundamental 
measure theory (FMT) for mixtures of general convex hard bodies |l4| . FMT has proved to account accurately for both 
the structure and the thermodynamics of inhomogeneous hard-sphere mixtures |l5| ] as compared with simulations. 
Recently, a FMT for a mixture of spheres and needles of vanishing thickness was proposed Our extension to 

include a nonsphcrical particle with finite volume (the rod) in the theory requires the deconvolution of the Mayer-/ 
function of the rod-sphere interaction into a set of orientation dependent weight functions |l4|] . The details of this 
technically involved deconvolution will be explained elsewhere ]l7| . 

We obtain the depletion potential in two steps. First we calculate the unperturbed density profile p s {z) of the 
spheres close to the wall. This is the equilibrium profile in the limit p r — > 0, as required in Eq.(l]). The functional 
reduces, in this case, to the well-tested FMT for spheres. In the second step we use p s (z) to evaluate c^\z, 9), using 
the deconvoluted Mayer-/ function. From Eq.([l]) we then obtain the depletion potential W(z, 9) for all z and 9 (from 
the profile p s (z) alone!). 

We present results for L/a = 10, a/a s = 1, and a packing fraction r/ s = (ir/6)a^p s = 0.2239 for the spheres. 
Figure || displays the depletion potential (3W{z 7 9) as function of z and 9. The first observation is that the potential 
is, unlike the AOa, not monotonic; the hard-sphere correlations generate a sequence of potential barriers and wells. 
The dashed and dotted curves in Fig. |^ denote the positions of the minima and maxima of the potential. We note 
that for small angles 9 <C 1 the shape of the rod-wall depletion potential W(z, 9) coincides almost perfectly with the 
depletion potential W ws (z) between a sphere of diameter a and the wall, i.e., W(z,6 « 1) ~ W ws (z — Lcos|0|/2). 
For these small angles the length of the rod L is rather unimportant for details of the depletion potential, as we have 
verified for various values of L/a and a/a s . For large angles, however, the whole geometry of the rod, i.e., L/a and 
a/a s , is relevant, and the depth of the depletion potential at contact becomes more negative as both L/a and a/a s 
are increased. In addition the contact value of the depletion potential can be further decreased by increasing rj s . 

Although certain general trends of the influence of the geometry on the shape of the depletion potential can already 
be roughly understood within the simple AOa [Q, the correlations in the hard-sphere fluid are very important and 
can lead to a quantitatively and even qualitatively different behavior. In Fig. ^ we illustrate this by comparing the 
depletion potential at contact as obtained within our DFT approach (W(z m i n (9), 9), full line) with the corresponding 
quantity within the AOa (W AOa (z m i n {9), 9), dotted line). While for small and large angles the contact value of the 
full depletion potential is more negative than within the AOa, for intermediate values of 9 packing effects can shift 
the actual contact value above its corresponding AOa- value. If the rod is almost parallel to the wall, i.e., for 9 < ir/2 
the deviation of the contact value calculated within the AOa from the prediction of the full DFT is significant and 
much more pronounced than in the well-studied wall-sphere or sphere-sphere geometry. The inset of Fig. |3j shows that 
a rod in contact with the wall can exhibit a metastable orientation for a certain angle < 9q < ir/2 while the AOa 
predicts that this configuration is unstable. Although for the parameters under consideration here the height of the 
potential barrier is rather small (0.007/cbT), we note that its value increases for smaller values of a/a s . 

As a consequence of the dependence of the rod-wall depletion potential on the orientation 9 for a given (fixed) 
distance z of the center of the rod from the wall an entropic torque acts on the rod which drives it into an orientation 
with minimal depletion potential. This torque can be estimated by replacing the spherocylinder by a dumbbell 
composed of two spheres with diameter a connected by an infinitesimally thin but rigid wire of length L. Since the 
sphere- wall depletion force f S w{z) acts on the two spheres located at z± = z ± (L/2) cos 9 the corresponding torque 
with respect to the center of mass is given by 

M db (z,9) = ^sme[f sw (z-)-f sw (z+)]- (2) 

It turns out that at low packing fractions of the hard-sphere solvent, or at sufficiently large separations of the rod 
from the wall Eq. (g) yields a semi-quantitative expression for the torque acting on a spherocylinder. 

From the full rod- wall depletion potential W(z, 6) one can obtain the torque by rotating the rod by an infinitesimal 
angle d9 around an axis through the center of the rod in a direction characterized by the unit vector rig normal to 
the symmetry plane shown in Fig. |l|. The corresponding change in the depiction potential is dW(z, 9) which can 
be written as dW — — ^ fi ■ dri, i.e., as a sum of forces fi acting on the rod at positions from the center with 
dr.i = d9{ri x ng). It follows that the torque M(z, 9) = M{z, 9)i\g with respect to the center of mass is given by 
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(3) 



The symmetry of the problem leads to M(z, 9) = for 9 = and 9 = ir/2. A positive value of the torque acts on the 
rod as to increase the angle 9 (rotating it parallel to the wall) while a negative value of M leads to a decrease of 9 
(rotating it normal to the wall) . Some typical examples for the torque as a function of 9 for various values of z are 
shown in Fig. |]. Lines and symbols denote results from DFT and molecular dynamics (MD) simulations, respectively. 
The DFT predictions are in excellent agreement with the simulations. The details of the simulations will be presented 
elsewhere fl7| . 

For larger distances of the rod from the wall W(z, 9) exhibits local minima and maxima in addition to the global 
minimum with corresponding zeros of the torque (see Fig. |J for z = 5.5ct s and z = 4.5cr s ). A minimum (maximum) of 
W(z, 9) leads to a zero of M(z, 9) with a negative (positive) gradient in 9. Furthermore, the torque can exhibit a cusp 
(see Fig. [|) if the minimal distance between the wall and the rod leaves space for precisely one solvent sphere. This can 
be realized for distances a s +a/2 < z < a s + (L+a)/2 and a corresponding orientation 9 cusp = arccos((2(z— a s )~ a)/L). 

As the rod moves closer to the wall the modulus of the torque increases and for small separations from the wall (see, 
e.g., z = 1.5er s in Fig. ^) the torque vanishes only for 9 = ir/2 which is, however, only a metastable configuration. For 
such small and fixed values of z the entropic torque rotates the rod towards configurations with smaller angles until 
the rod is in contact with the wall and has to stop its rotation. Thus the rod reaches its most favorable configuration 
of lying parallel in contact with the wall not by approaching the wall in a parallel configuration but by touching the 
wall first at one end and then by decreasing the distance of its center from the wall. The modulus of the maximum 
of the entropic torque for the system considered here is of the order of about 20fcsT rad -1 , which corresponds to 
roughly 10~ 20 J rad -1 at room temperature. The strength of the torque increases for larger values of L/a, a/a s , or 
rj s . For example, for L/a = 20, a/a s = 2, and % = 0.2239, the modulus of the maximum of the entropic torque 
reaches a value of about 38.7fcsT rad -1 . We note that the maximum torque acts at small values of z for which the 
AOa also would predict the existence of a torque; however, its magnitude is largely enhanced through the correlations 
in the solvent. 

With the results presented so far we can comment on some aspects of the path in the z-9 plane a rod would take 
upon approaching the wall from the bulk. Only if the rod comes sufficiently close to the wall to be subject to the 
oscillations of the number density of the hard-sphere fluid as a function of z, the entropic force and torque will act on 
it. The closer the rod gets to the wall the higher are the potential barriers between the minima (see Fig. ||) the rod 
has to overcome by thermal motion in order to move still closer. The potential barriers are, however, only moderate 
and easy to overcome for small angles and increase for angles 9 > 60 degrees so that it is not optimal to rotate the rod 
until it is parallel to the wall since then the potential barrier for further approach is largest and for the parameters 
considered here more than 2.3kgT. On the other hand along the line of the second minimum - see the dashed lines 
in Fig. ||- (3W has a similar shape as (3W(z m i n {9),9) [see Fig. [|, so that the entropic torque will rotate the rod into 
the local minimum at 9 = ir/2 if 9 is already sufficiently large. Once the rod is in contact with the wall the torque 
will rotate it towards the pronounced global minimum of /3W(z, 9). 

We have presented quantitative predictions of the entropic torque acting on a hard sphcrocylinder close to a 
hard planar wall in a solvent of small hard spheres. Our DFT predictions are in excellent agreement with our MD 
simulations. We find that the the depletion effect leads to a significant entropic torque and it is tempting to speculate 
that this entropic torque can play an important role in understanding the key-lock principle. 
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FIG. 2. Depletion potential (3W(z, 9) for the spherocylinder shown in Fig. y for L/a = 10, a/a s = 1, and a packing fraction 
rj s — p s 7ra"s/6 = 0.2239 of the hard-sphere solvent. 8 is measured in degrees. The dashed and dotted lines represent the 
positions of local minima and maxima, respectively. 
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FIG. 3. The contact value [3W (zmin{8) , 8) as a function of 8 (in degrees) as calculated within DFT (full line) and within 
the AOa (dotted line) for the system described in Fig. Q The inset shows that a rod in contact with the wall can exhibit a 
metastable orientation for 8 = 8q ~ 8 degrees with a barrier height /3AW = 0.007. 
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FIG. 4. The torque (3M(z, 8) as function of 8 for various values of z. For reasons of clarity the curves for z = 4.5, 3.4, . . . 
are shifted downward by —5, —10, . . ., respectively. The horizontal lines M — are indicated partially. The parameters of the 
system are the same as in Fig. |^. Because of the hard wall at z — the torque is defined only for 8 > arccos(2z/L — cr/L) if 
z < (a + L)/2. M > (< 0) corresponds to a force which tends to align the rod parallel (normal) to the wall. The lines denote 
our DFT results which are in excellent agreement with our simulations (symbols). The error bars of the simulations are of the 
order of the symbol size. 
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